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iy~j Abstract: Surface waves in classical fluids experience a rich array of 

<0 black/white hole horizon effects. The dispersion relation depends on the 

characteristics of the fluid as well as on the fluid depth and the wavelength 

regime. We focus on the shallow-water regime, and discuss the experimental 

f— V proof that the circular hydraulic jump marks the transition between a su- 

CN percritical and a subcritical flow regime. This finally confirms a theoretical 

'Tt conjecture formulated by Lord Rayleigh nearly 100 years ago. It also con- 

>^ firms that the circular jump corresponds to the spontaneous formation of 

S^ a hydrodynamic white hole, with interesting characteristics from the point 

H of view of analogue gravity. We study the dispersive regime, mention some 

lessons about the trans-Planckian issue and describe possible directions for 

future work. 
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1.1 Introduction 

Surface waves in classical fluids provide a natural and rich class of black/white 
hole analogues. A familiar example of a hydrodynamic black hole analogue 
is the draining bathtub. From an experimental point of view, this has the 
obvious complication that the interesting physics literally disappears down 
the drain, and is unavoidably mixed with effects such as vorticity and tur- 
bulence which complicate the analysis. But many properties of black holes 
are directly related to those of white holes through a simple time-inversion, 
at least in the non-dispersive and non-dissipative limit. Experimental ef- 
forts have therefore so far concentrated mainly on white hole analogues. 
Two familiar examples are the blocking of sea waves at a river mouth and 
the approximately circular jump created by opening the tap in a kitchen 
sink. These two types of white hole analogues can be reproduced in con- 
trolled laboratory settings in order to study the associated horizon effects 
and their possible lessons for relativity (and vice versa: lessons from relativ- 
ity for fluid mechanics). The river- mouth example corresponds to deep water 
waves propagating against a counter-current in a wave channel, while the 
kitchen-sink example corresponds to shallow waves on a circular hydraulic 
jump. 

The general dispersion relation for capillary-gravity surface waves prop- 
agating against a counter-current of velocity U (see, e.g., [1] and references 
therein) is 

{uj - V.kf = (gk + ^kA tanh{kh), (1.1) 

with g the gravitational constant, 7 the surface tension, p the density, h the 
fluid depth, and as usual uj and k are the frequency and wavenumber, respec- 
tively. From now on we will neglect any tangential effects and consider only 
surface waves propagating linearly (or radially) against the background flow. 
The "relativistic" regime is obtained in the shallow- water limit (kh ^ 1) of 
gravity- wave propagation (7 = 0): (w — Uk)'^ ~ ghk'^, with the relativistic 
"invariant" speed c = -v/^/i [2j . 

It is also instructive to study deviations from this relativistic regime, 
where the propagation velocity of the surface waves is no longer invariant. 
From now on, we will therefore refer to a "horizon" in the generalized sense 
of the point (or region) where the group velocity Cg = duj/dk of the surface 
waves vanishes. Indeed, such deviations can provide useful clues with regard 
to quantum gravity phenomenology and the robustness of black hole physics 
(for instance. Hawking radiation) with respect to possible non-relativistic 
features of the trans-Planckian physics. Besides Hawking radiation, other 
applications include the study of quasinormal modes (QNMs): the decay or 
energy dissipation modes of a perturbed black-hole spacetime. Such QNMs 
can be studied analytically and numerically, for example for surface waves 
(or sound waves) in the draining bathtub geometry [3i|4j. "Trans-Planckian" 



corrections to this QNM spectrum could then be similar to those obtained 
for acoustic black holes in Bose- Einstein condensates [5], since the high- 
wavenumber dispersion relation of surface waves has a superluminal /c^- 
dominance similar to the BEC dispersion relation 

(a; - Ukf = cl {k^ + \ek'^ , (1.2) 

where Cg is the speed of sound and ^ = h/{mcs) represents the healing length 
of the condensate. 

We will now focus on a typical example of a shallow-water geometry, 
namely the circular hydraulic jump. 



1.2 The circular hydraulic jump 

When a vertical fluid jet impacts on a horizontal plate with a sufficient flow 
rate, it will form a thin layer near the impact zone, which expands radially 



and at a certain distance abruptly increases in thickness, see Figs. 1.1 and 
1.2[ This phenomenon is known as the circular hydraulic jump, and was 
already described by Leonardo Da Vinci. Its first modern description is 
due to Lord Rayleigh [6], who developed a momentum-balance theory to 
describe it, and hypothesized that its appearance is due to a deceleration 
from a supercritical flow in the inner region to a subcritical flow outside. In 
other words, in the inner region, the radial fluid velocity at the surface v^ 



(which plays the role of the counter-current U in the dispersion relation ( 1.1 ) 
above) is larger than the propagation speed of the surface waves c, so that 
these surface ripples can only propagate downstream. In the outer region, 
v^ < c and hence the ripples can propagate in both directions. The transition 
would be characterised precisely by the jump, analogous to a shock wave. 
The circular jump is therefore the stationary wave solution or zero mode 
(zero- frequency solution of the dispersion relation) in the case of a circularly 
propagating, decelerating fluid flow. 

The essential point in the context of analogue gravity is similar to the 
case of acoustic black holes [3 [8] and other examples of analogue gravity [9] : 
the propagation of these surface waves obeys a generalised d'Alembertian 
equation in which the intervening curved-spacetime metric is the 2-|-l-di- 
mensional Painleve-Gullstrand-Lemaitre (PGL) form of the Schwarzschild 
metric [lOJ. Indeed, the line element for the circular jump is 

ds^ = 1 [[c^ - {vlf]dt^ + 2vldt dr - dr^ - r^dc/)^] , (1.3) 

where the surface wave propagation speed c plays the role of the speed of 
light in gravity, while the radial surface flow velocity v^ corresponds to the 




Figure 1.1: The circular hydraulic jump. 




Figure 1.2: Basic setup of the circular hydraulic. A liquid is pumped through 
a nozzle and the fluid jet impacts vertically onto a horizontal plate. At a 
distance Rj, a sudden "jump" occurs in the fluid height. 



local velocity of an observer freely falling from infinity in the case of gravity. 
The PGL metric 

--ft'' 5) e-^) 

with / the unit matrix, transforms into the Schwarzschild form through 
the coordinate transformation di = dt + drv^/[c^ — {v^)"^]. The condition 
gtt = for a horizon in the non- rotating case becomes simply c = v^. It 
is in this precise sense that the circular jump is believed to constitute the 
hydrodynamical analogue of a white hole |1H [T2] . 

Ever since Watson improved Rayleigh's account by incorporating the 
effects of viscosity p!3], there has been a continuous experimental and theo- 
retical interest in the circular jump, with a remarkable burst of work in the 
last decade. Curiously, though, an experimental proof of Rayleigh's original 
conjecture that the jump corresponds to a supercritical-to-subcritical transi- 
tion (and hence, in the language of analogue gravity, that the circular jump 
constitutes a white hole) was only provided recently through the Mach cone 
experiment |14j that we describe next. 



1.3 Mach cone experiment 

The question of whether and where the surface waves are blocked by the 
fluid flow could be answered by comparing the radial fluid velocity v], at the 
surface and the propagation velocity c of the surface waves. However, given 
that we are not interested in the precise values of each quantity separately, 
but only in the ratio v^/c, we have opted for the following alternative. Our 
method is based on the Mach cone well known in the case of sound waves. 



see Fig. L3 There, an object propagating at a speed V above the speed 
of sound Cs leaves behind an observable cone, the Mach cone. The cone is 
formed by the envelope of the subsequent wavefronts emitted by the object, 
which partially escape from each other. The half-angle 9 of the cone obeys 



sin 6* = ct/\V\t = c/\V\ = 1/M, with M the Mach number, see Fig. [O 
For an object propagating at a speed V < Cg-, the subsequent wavefronts 
remain inside the previous ones and no Mach cone is formed. Note that 
the appearance of a Mach cone can be interpreted as a consequence of the 
Cerenkov threshold condition expressed as w — k.V = 0, or in other words 
uj — k.V cos ip = 0. Since the Cerenkov angle "0 and the Mach angle 9 are 
complementary, i.e. ip = tt/2 — 9, the same trigonometric rule sm9 = c/\V\ 
again follows immediately. 

Exactly the same argument can be applied to the case of an object 
standing still at the surface of a fluid flow. By comparing the fluid flow 
velocity v'^ to the propagation speed of surface waves c, the following cases 
can occur: 



• 



Supercritical region: v^ > c =► ^ G [O, |] 



Subcritical region: v^ < c 



6 complex; Mach cone disappears 



• Horizon: v^ = c 



Subsonic |v| < c 




Supersonic |v| > c 




Mach Line 



Figure 1.3: Wavefronts emitted in subsonic/subcritical regime (left) and 
formation of Mach cone in supersonic/supercritical regime (right). 




Figure 1.4: Circular hydraulic jump and Mach cone created by inserting a 
sharp object in the fluid flow. 

Our experiment to demonstrate the presence of a hydrodynamic horizon 
is described in detail in [T3]. Essentially, we have pumped silicon oil through 
a steel nozzle onto a horizontal PVC plate. A needle was placed such as 
to penetrate the flow surface at varying distances r from the centre of the 
circular jump, where the oil jet impacts on the PVC plate. For each position 
of the needle, we haved photographed the setup with a high-speed camera 



and measured the corresponding Mach angle 9, see Fig. 1.5 The resulting 



angles and the derived ratio v^/c are shown in Fig. 1.6 Inside the jet impact 
zone (r < a, where a is the nozzle radius, see also Fig. 1.2 above), we expect 



<C c followed by a steep increase for r > a until a certain value u 



s,max 

'r , 



since the fluid impacts vertically before being converted into a radial flow. 
The field of vision of our experimental setup starts near this maximum, see 



Fig. 1.6, corresponding to a Mach angle 9 of roughly vr/lO. From there, 
smoothly increases to about 7r/4 at approximately 3/4 of the jump radius, 
and then rapidly opens up to reach exactly 7r/2, and hence v^/c = 1, near 
the ridge of the jump. 
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Figure 1.5: Measurements of the Mach angle 6. A needle is placed inside 
the flow at varying distances from the centre of the jump, (a) Mach cone 
near the centre of the jump, (b) Mach cone near the edge of the jump, (c) 
The Mach cone disappears just outside the jump. 

These results provide a clear proof that the circular hydraulic jump con- 
stitutes a two-dimensional hydrodynamic white hole: surface waves travel- 
ling at a velocity c towards the jump from the exterior are trapped outside in 
precisely the same sense as light is trapped inside a gravitational black hole. 
The corresponding white-hole horizon is situated at the radius of the jump 
itself. Curiously though, in our experiments, the critical point v^/c = 1 is 
actually reached in a very smooth way. The ratio v^/c decreases mainly 
far from the jump, well inside the inner region. This is in striking contrast 
with the standard theoretical models in fluid mechanics, which describe the 
circular jump as a shock wave and therefore prescribe that the critical point 
itself should lie within a sharp (and in models without viscosity: discontin- 
uous) transition from a supercritical to a subcritical regime, see e.g. |15j. It 
is not clear whether this smooth transition to the critical point is a genuine 
property of the jump itself, or a consequence of the perturbation of the flow 



100 



30 - 



60 



40 



20 



. M*-***-^ 



I 
I 

i 

M 

\ 



y 






10 



20 

r (nim| 



> 



30 p 40 



Figure 1.6: Mach angle 9 (red triangles) and ratio vf./c (blue circles) as a 
function of the distance r from the centre of the jump. The dashed vertical 
lines represent the nozzle radius a and the jump radius Rj. Experimental 
parameters: see [TIJ. 



pattern due to the insertion of the needle, which somehow smoothens out 
the shock wave (or shifts its position). But the main result stands clear: the 
circular jump is a white hole. 

From the point of view of analogue gravity, the following features of this 
hydrodynamic white hole are particularly striking, especially in comparison 
with other current or planned experiments, e.g. in optics [16] or Bose- 
Einstein condensates |17j: 

• The white hole here is created "spontaneously". One only needs to 
arrange for a sufficient fluid flow rate, but no extraordinary engineering 
is required to obtain a long-lived horizon. In fact, it suffices in principle 
to open the tap in a kitchen sink to observe such a white hole. All 
other ingredients in the experimental setup (the choice of silicon oil, 
the precision of the pump etc) serve merely to make the experiment 
cleaner and free of perturbations and other undesired effects, but do 
not affect the essential point which is the spontaneous creation of a 
(hydrodynamic) white hole. This makes hydrodynamic white holes 
a rather exceptional example within the condensed-matter analogues, 
as in most analogue models there is a strong resistance against the 
creation of a horizon, see e.g. |18j . which is usually accompanied by 
the appearance of several types of instabilities which tend to destroy 



the horizon. 

• The white hole can be observed with the naked eye: the location 
of the white hole horizon is precisely where the fluid undergoes the 
characteristic jump. This can be verified in the way described above, 
in principle simply by sticking a pencil in a kitchen sink. 

These two elements mean that the circular jump provides a nice example of 
a simple "do-it-yourself black-hole experiment" which can be explained in 
terms accessible to any person with a basic knowledge of trigonometry. 

Let us now go back to the dispersion relation for surface waves and 
examine in which dispersive regime the circular jump typically lies. 



1.4 Dispersion relation 

Developing Eq. ( |1.1| ) for the case kh <^ 1 and truncating at 0{k'^), one 
obtains 



(w - Uky ^ ghk^ + ( V ) ^^ ^^^^ ^"^'^^ 

= c^k' + c^fp^-^^k'^ + Oik^), (1.6) 

where Ic = \/"f/gp is the capillary length {Ic = 1.49mm for the silicon oil 
in our experiments). The shallow- water regime is relativistic at low values 
of k (the associated "relativistic speed" c is c = ^/gh, i.e. the surface wave 
velocity in the low-A; or gravity limit where capillarity is negligible). For in- 
termediate fe's, the dispersion can either be "subluminal" or "superluminal" 
(i.e., the group velocity Cg = -J^ decreases or increases with k, respectively). 



see Fig. 1.7, with a critical transition depth /itrans = V^lc (see also [T9]). 
corresponding to an inflection point cPcg/dk'^ = at A; = 0. At higher fe's 



(once /c/i ~ 1), the development (1.5) can no longer be trusted, and the dis- 
persion becomes superluminal, irrespective of the value of h, asymptoting 
to Cg DC vfc- 

We are here interested in the intermediate regime, i.e. in the first non- 
relativistic corrections. For our silicon oil, /itrans ~ 2.6mm, while the typical 
heights of the fluid inside the jump are much smaller (typically ~ 1mm). 
The same relation h < \/3lc holds true for experiments with other fluids 
reported in the literature (at least in the so-called "type I" circular jump, 
which avoids reversal currents and is therefore the most interesting from 
the point of view of analogue gravity, see [13] )• This confirms that the 
circular jump is superluminal {cg increasing with k) for realistic experimental 
parameters. A direct observational test is the following. When introducing 
the needle in a circular jump inside a water flow, tiny capillary waves are 



observed upstream from the needle. In other words, these capillary waves 
manage to propagate some distance against the background flow, thereby 
confirming the superluminal character. Note that in water, /itrans ~ 4.7mm 
whereas typical fluid depths inside the jump are < 0.5mm |20j . In silicon oil, 
such capillary waves upstream from the needle are not observed, indicating 
that the superluminal correction is too small to compensate for the strong 
dissipation due to viscosity. 



fj (rad/s) 
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Figure 1.7: Group velocities (at U = 0) for silicon oil (viscosity i^ = 20cS ss 
SOi^water; surfacc tension 7 = 0.0206N/m ss ^7watcr; density p = 950kg/m^). 
and different values of the fluid height h [from bottom to top at low u]: 
h = 1mm (red dashed line); 2.5mm (green dashed-dotted line); 3mm (blue 
dotted line); 5mm (full brown line). Typical fluid heights in the experiments 
reported in |T3j are ~ 1mm, corresponding to the red dashed line and hence 
a purely superluminal correction to the dispersion relation. 

In such a superluminal regime, sufficiently high-frequency modes can 
penetrate across the horizon in the classically prohibited sense. The interior 
of the black/ white hole is then no longer causally separated from the outside. 
Such superluminal horizon-crossing effects are considered in several scenarios 
for quantum gravity phenomenology. They are especially interesting in the 
sense that they put the robustness of several aspects of black hole physics 
(and in particular, Hawking radiation) to the test [2l]. Most of these are 
ultimately related to vacuum selection issues [HI |23] , which are absent in 
the case of the circular jump: There is a clear "preferred" reference frame, 
the laboratory reference frame, and hence a clear and well-defined vacuum 
state (the one in which the fluid stands still with respect to the laboratory). 
Nevertheless, some aspects of the robustness of the thermal spectrum of 
Hawking radiation with respect to trans-Planckian physics could in principle 
be tested in the circular hydraulic jump, see the discussion in [14J . 

The main barrier to such studies is the fact that these penetrating modes 
are blueshifted by the countercurrent, and therefore strongly damped due 
to viscosity, making it extremely hard to measure them in a reliable way. 
Superfluids could provide an interesting way out of this problem. In fact, 
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an earlier experiment with a circular jump in ^He |24j has shown capillary 
ripples near the jump. Below the superfluid transition temperature, station- 
ary capillary ripples also appeared well within the inner region. There also 
exists a concrete proposal for an experiment with ripplons at the interface 
between the superfluid ^He-A and ^He-B phases in a rotating cryostat [25], 
related to an earlier experiment reported in [26]. The essential idea is to 
take advantage of the fact that the critical rotation velocity Oc for the ap- 
pearance of vortices, which allow the superfluid component to co-rotate with 
the normal component, is lower for the A-phase than for the B-phase. A 
cylindrical sample with superfluid ^He is set into rotation, and separated 
into the A- and B-phase by the gradient of a magnetic field. A differential 
rotation 0, between both superfluid phases is achieved at rotation velocities 
ftc,A < ^ < ^c,B- By modulating the propagation velocity c of the ripplons 
through a variation of the thickness of the B-layer channel, the interface 
flow velocity v = rCl can be made to locally exceed c (typically of the order 
of 1 mm/s), thereby creating an ergoregion where v > c. For sufficiently 
thin B-channel layers, the relativistic limit corresponding to the shallow- 
water regime described earlier is recovered, and the boundaries w = c of the 
ergoregion form a pair of relativistic black/white hole horizons. 

Such an experiment could be interesting not only to measure the radi- 
ation associated with these horizons, but also because it can be extended 
to simulate a black hole singularity. Indeed, at an even higher critical ve- 
locity, corresponding to the proper Kelvin-Helmholtz threshold [25] , the gu 
component of the effective metric becomes infinite, indicative of an effective 
black hole singularity. 

The question of dissipation is also related to the trans-Planckian issue, 
which we discuss next. 

1.5 Trans-Planckian problem 

In analogue gravity systems, the infinite blue-shifting associated to strictly 
relativistic horizons is avoided because of dispersive corrections at short 
wavelengths. One of the cornerstones of the analogue gravity programme 
has therefore from the start been that perhaps the same could occur in real 
gravity. 

Historically, the focus has been mainly on the study of subluminal ( "nor- 
mal") dispersion. In that case, the trans-Planckian problem is avoided be- 
cause the incident waves are reflected just before reaching the white hori- 
zon, since their group velocity decreases as they blueshift. For a black hole 
horizon, the argument should be time-reversed: late-time outgoing modes 
originate from finite-A: modes which start off outside the horizon and are 
redshifted as they are dragged towards the black hole. Due to this (finite) 
redshifting, their group velocity increases and they manage to escape away 
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from the black hole. 

Superluminal dispersion relations, on the other hand, solve the trans- 
Planckian problem by allowing blue-shifted waves to penetrate across the 
(low-fc) horizon, rather than endlessly blueshifting against it, as in the purely 
relativistic case. However, if one imagines a flow with a counter-current 
which increases without bounds, then the superluminal waves will also be 
blueshifted without bounds as they propagate ever further across the hori- 
zon, against an ever faster counter-current. Again, for a black hole, the 
argument should be time-reversed, and outgoing superluminal modes seem 
to ultimately originate from the singularity. 

An additional characteristic scale of the system is needed to avoid this 
scenario. In principle this could be some scale which puts a limit on the 
acceleration of the fluid. Once the counter-current reaches its maximum 
value, the superblueshifted capillary waves would then be allowed to freely 
propagate against the background flow. In practice, however, in the case of 
the circular jump in a classical liquid, the characteristic scale which avoids 
the endless blueshift is associated to dissipation: the blueshifted waves are 
strongly damped due to viscosity, typically very close to the horizon (just 
inside the jump). The trans-Planckian problem is thus solved through a 
combination of dispersion and dissipation. The analogue of the Planck scale 
(the intermolecular distance) is never reached and the fluid continuity ap- 
proximation is never violated. 

It is an open question how the trans-Planckian problem is solved in 
real gravity, if at all. No clues of modified dispersion or other forms of 
Lorentz violation have been found so far, and there exist strong bounds on 
the possibility of Lorentz violations at scales at least up to lO^E'pianck, see 
e.g. |271 128| . The question is then whether this means that trans-Planckian 
frequencies are really reached in the physics associated to black holes, imply- 
ing that some theory of quantum gravity comes into play, or whether there 
exist other (non-dispersive) mechanisms to avoid such a trans-Planckian 
scenario. The latter was argued, for example, in [29]. Assuming that black 
holes in the strict sense of general relativity (with long- lasting trapping hori- 
zons and a central singularity) do not form in astrophysics, but are replaced 
by "black stars", objects which hover just outside their own gravitational 
radius (maintained, for example, by vacuum polarization pressure [30j), it 
was argued that Hawking radiation could be produced by the internal pul- 
sations of such objects without any appeal to trans-Planckian frequencies 
nor the need to invoke modifications of the dispersion relation. 

1.6 Outlook 

The proof that the circular jump constitutes a hydrodynamic white hole 
confirms Rayleigh's century-old conjecture that the circular jump marks the 
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transition between a supercritical and a subcritical regime. The proof itself 
was based on well-known physics (the Mach cone) and used a straightfor- 
ward methodology. This relative simplicity means that it can in principle 
be repeated in a kitchen sink. But this should not deceive anyone: further 
advance in the study of the circular jump as a hydrodynamic white hole is 
likely to be much more complicated and involve state-of-the-art technology. 
Detecting the radiation associated with this hydrodynamic horizon would 
require either precision-measurement of the blueshifted waves that appear 
inside the jump, or the use of viscous-free fluids such as liquid Helium near 
zero temperature. Whether such approaches will allow results similar to 
those obtained in the wave-tank experiments j31| !32| is an interesting ques- 
tion. The complication lies mainly in the dissipation of the blueshifted waves 
near the horizon. The circular jump has a double disadvantage. The first 
one is a matter of scales: it is easier to work with long wavelengths, which 
are less prone to dissipation, in a linear water-tank than in a circular set- 
up. Second, and associated: the superluminal character of the dispersion 
relation in the case of the circular jump means that the Hawking process 
involves the physics inside the horizon, where there is additional blueshifting 
due to the accelerating counter-current, whereas in the subluminal regime 
of the water-tank, all the processes related to Hawking radiation take place 
outside the horizon. 

The question of whether it is possible at all to reproduce the hydro- 
dynamic analogue of Hawking radiation in a situation which involves the 
physics inside the horizon, would also precisely make the circular jump an 
interesting test case. Several "trans-Planckian" issues involved and compet- 
ing effects were discussed in [Uj. 

Whatever the outcome of such further work, it should be clear that the 
interest of cross-fertilization between fluid mechanics and gravity, contrarily 
to a relativistic horizon, is two-directional. 
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